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Abstract 

The paper deals with an approach to the model-independent searching 
for the Z' gauge boson as a virtual state in scattering processes. The 
relations between the Z' couplings to fermions covering a wide class of 
models beyond the standard model are found and used. They reduce in 
an essential way the number of parameters to be fitted in experiments. 
Special observables which uniquely pick out the Z' at energies of LEP and 
ILC colliders in different leptonic processes are introduced and the data of 
LEP experiments are analyzed. The Z' couplings to leptons and quarks 
are estimated at 95% confidence level. At this level, the LEP data are 
compatible with the existence of the Z' with the mass m z > ~ 1 — f .2 TeV. 
These estimates may serve as a guide for experiments at the Tevatron 
and/or LHC. A comparison with other approaches and results is given. 



1 Introduction 

The precision test of the standard model (SM) at the LEP gave a possibility 
not only to determine all the parameters and particle masses at the level of 
radiative corrections but also afforded an opportunity for searching for signals 
of new heavy particles beyond the energy scale of it. On the base of the LEP2 
experiments the low bounds on parameters of various models extending the 
SM have been estimated and the scale of new physics was obtained [TJ [51 [31 
[H [5] . Although no new particles were discovered, a general believe is that the 
energy scale of new physics to be of order 1 TeV, that may serve as a guide for 
experiments at the Tevatron and LHC. In this situation, any information about 
new heavy particles obtained on the base of the present day data is desirable 
and important. 

Numerous extended models include the Z' gauge boson - massive neutral 
vector particle associated with the extra U{\) subgroup of an underlying group. 
Searching for this particle as a virtual state is widely discussed in the literature 
(see Ref. [SI [7] for review) . In the content of searching for Z' at the LHC and the 
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ILC an essential information and prospects for future investigations are given 
in lectures [8]. Such aspects as the mass of Z', couplings to the SM particles, 
Z-Z' mixing and its influence in various processes and particles parameters, 
distinctions between different models are discussed in details. We shall turn 
to these papers in what follows. As concerned a searching for Z' in the LEP 
experiments and the experiments at Tevatron [5], it was carried out mainly in 
a model-dependent way. Some popular models has been investigated and low 
bounds on the mass mz< were estimated (see Refs. [TJ [21 [21 SJ [S] , the recent 
results in Refs. PHHH]). As it is occurred, the low masses are varying in a wide 
energy interval 400-1800 GeV dependently on a specific model. These bounds 
are a little bit different in the LEP and Tevatron experiments. In this situation 
a model-independent analysis is of interest. 

In the papers |T2l [T3j [XSJ [16] of the present authors a new approach for 
the model-independent search for Z'-boson was proposed. In contrast to other 
model-independent searches, it gives a possibility to pick out uniquely Z' virtual 
state proper for a wide class of models (listed below) beyond the SM. Our con- 
sideration is based on two constituents: 1) The relations between Z' couplings 
motivated by renormalizability of an unknown theory beyond the SM. Due to 
these relations, a number of unknown Z' parameters entering the amplitudes 
of different scattering processes considerably decreases. 2) When these rela- 
tions are accounted for, some kinematics properties of the amplitude become 
uniquely correlated with this virtual state and the Z' signals exhibit themselves. 
The corresponding observables have also been introduced and applied to ana- 
lyze the LEP2 experiment data. Comparing the mean values of the observables 
with the necessary specific values, one could arrive at a conclusion about the 
Z 1 existence. The confidence level (CL) of these values has been estimated and 
adduced in addition. Without taking into consideration the relations between 
coupling the determination of Z'-boson requires a supplementary specification 
due to a larger number of different couplings contributing to the observables. 

In Refs. [T^l H51 HH] the one-parametric observables were introduced and 
the signals (hints in fact) of the Z' have been determined at the la CL in the 
e + e~ — > /i + /J~ process, and at the 2a CL in the Bhabha process. The Z' mass 
was estimated to be 1-1.2 TeV. An increase in statistics could make these signals 
more pronounced. In Ref. |17] the updated results of the one-parameter fit and 
the complete many-parametric fit of the LEP2 data were performed with the 
goal to estimate a possible signal of the Z'-boson with accounting for the final 
data of the LEP collaborations DELPHI and OPAL [2 [3j IU E] . Usually, in a 
many-parametric fit the uncertainty of the result increases drastically because of 
extra parameters. On the contrary, in our approach due to the relations between 
Z' couplings there are only 2-3 independent parameters for the investigated 
leptonic scattering processes. As it was showed in Ref. [T7], an inevitable 
increase of confidence areas in the many-parametric space was compensated 
due to accounting for all accessible experimental information. Therefore, the 
uncertainty of the many-parametric fit was estimated as the comparable with 
previous one-parametric fits in Refs. |151 116[ . In this approach the combined 
data fit for all lepton processes is also possible. Note that the hints for the Z' 
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have been determined in all the processes considered that increases the reliability 
of the signal. These results may serve as a good input into the LHC and future 
ILC experiments and used in various aspects. To underline the importance of 
them we mention that there are many tools at the LHC for the identification 
of Z' . But many of them are only applicable if Z 1 is relatively light. The 
knowledge of the Z' couplings to SM fermions also have important consequences. 
As concerns the notion "model-independent search" used below, it refers to a 
class of models containing Z 1 and inspired by the grand-unified field theories. It 
does not mean all possible ones. But if one determines the signal of this state, 
further specification of the underlying theory could follow. 

The paper is organized as follows. In sect. 2 we give a necessary information 
about the description of Z' at low energies and introduce the relations between 
the Z' couplings. In sect. 3 the cross sections and the observables to pick out 
uniquely the virtual Z' in the e + e~ — > [i + ,t + t~ processes are given. The 
fits of data are described and discussed. Then in sect. 4 the same is present for 
the Bhabha process e + e~ — > e + e~ . The one parametric and two parametric fits 
are discussed. In sect. 5 we discuss the role of the present model-independent 
analysis for the LHC experiments. The discussion and comparison with results 
of other approaches are given in sect. 6. 



2 The Abelian Z' boson at low energies 

Let us adduce a necessary information about the Abelian Z'-boson. This particle 
is predicted by a number of grand unification models. Among them the Eq and 
SO(10) based models [18] (for instance, LR, x~ and so on) are often discussed 
in the literature. In all the models, the Abelian Z'-boson is described by a low- 
energy U(l) gauge subgroup originated in some symmetry breaking pattern. 

At low energies, the -2T'-boson can manifest itself by means of the couplings 
to the SM fermions and scalars as a virtual intermediate state. Moreover, the Z- 
boson couplings are also modified due to a Z-Z' mixing. In principle, arbitrary 
effective Z' interactions to the SM fields could be considered at low energies. 
However, the couplings of non-renormalizable types have to be suppressed by 
heavy mass scales because of decoupling. Therefore, significant signals beyond 
the SM can be inspired by the couplings of renormalizable types. Such cou- 
plings can be derived by adding new [/(l)-terms to the electroweak covariant 
derivatives D cw in the Lagrangian [19l [20] (review, Ref. [HE]) 
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where summation over all the SM left-handed fermion doublets, leptons and 
quarks, f L = (/„)/,, {fd)h, and the right-handed singlets, f R = 
is understood. In these formulas, g,g',g are the charges associated with the 
SU(2)l, U(1)y, and the Z' gauge groups, respectively, er Q are the Pauli ma- 
trices, Qf denotes the charge of / in positron charge units, is the U(l)y 
hypercharge, and Yf L = — 1 for leptons and 1/3 for quarks. In general, genera- 
tors Yf L — dia,g(Yf u ,Yf d ) and Y^ = diag(Y0 ) i , ^,2) are diagonal 2x2 matrices. 
As for the scalar sector, the Lagrangian can be simply generalized for the case 
of the SM with two light Higgs doublets (THDM). 

The Lagrangian (fTJ) leads to the Z-Z' mixing. The Z-Z' mixing angle 9q is 
determined by the coupling Y$ as follows 

grinflWcos^^ i) + o(<), (3) 
V4na cm m l z , v \m%, J 

where Ow is the SM Weinberg angle, and a cm is the electromagnetic fine struc- 
ture constant. Although the mixing angle is a small quantity of order m^, , it 
contributes to the Z-boson exchange amplitude and cannot be neglected at the 
LEP energies. 

In what follows we will also use the Z' couplings to the vector and axial- 
vector fermion currents defined as 

<\i — ^ " f ^ • W 

The Lagrangian ([2]) leads to the following interactions between the fermions and 
the Z and Z' mass eigenstates: 

C zff = [{vfi + 7 5 a^) cos 9 + (v f + J 5 a f ) sin 9 ] f, 

Cz'ff = \z'jY[{vf +l b a f )cos9 -{vfi + j 5 a™) sin 9 ] f , (5) 

where / is an arbitrary SM fermion state; v^, ar^ are the SM couplings of 
the Z-boson. 

At low energies the Z' couplings enter the cross-section together with the 
inverse Z' mass, so it is convenient to introduce the dimensionless couplings 

m z _ m z , 
a f = —j= a /> v f = —j= v h ( 6 ) 

which can be constrained by experiments. 

The low energy parameters Y^, Ylj, Yrj must be fitted in experiments. In 
most investigations they were considered as independent ones. In a particular 
model, the couplings Y^, Ylj, Yrj take some specific values. In case when 
the model is unknown, these parameters remain potentially arbitrary numbers. 
However, this is not the case if one assumes that the underlying extended model 
is a renormalizable one. 
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In Refs. |T2] |T3] it was shown that these parameters are correlated. This cor- 
relations follows if the underlying unknown theory is a renormalizable one. The 
following conditions were assumed to derive the relations between Z' couplings: 

1. only one neutral vector boson exists at energy scale about 1-10 TeV, 

2. the Z' boson can be phenomenologically described by the effective La- 
grangian ©, ([J) at low energies, 

3. the Z' boson and other possible heavy particles are decoupled at consid- 
ered energies, and the theory beyond the Z' decoupling scale is either one- 
or two-Higgs-doublet standard model, 

4. the SM gauge group is a subgroup of possible extended gauge group of the 
underlying theory. So, the only origin of possible tree-level Z' interactions 
to the SM vector bosons is the Z-Z 1 mixing. 

Under these conditions, we have obtained the relations between phenomenolog- 
ical parameters of the effective Lagrangian ([2]), ([]}: 

%,x=%,i = %, Y LJ =Y LJ *, Y R j=Y LJ + 2Tz f Y <l> , (7) 

Here / and /* are the partners of the SU(2)l fermion doublet (I* — v\,v* = 
Z,g* = qd and q* d — q u ), T^f is the third component of weak isospin. They are 
key point for investigations reported below. 

Introducing the Z' couplings to the vector and axial-vector fermion currents 
(@J, the last formula in Eq. ([7]) yields 

v f -a f = v f , - a f * , a f = T 3f g%. (8) 

The couplings of the Abelian Z' to the axial-vector fermion current have a 
universal absolute value proportional to the Z 1 coupling to the scalar doublet. 
Then, the Z-Z' mixing angle ([3]) can be determined by the axial- vector coupling. 
As a result, the number of independent couplings is significantly reduced. 

We assume no new light particles. The relations could change essentially if 
the SM has to be modified at energies below the Z' mass. 

The derived relations are necessary but not exhaustive constraints on the 
Z' couplings. To derive exhaustive constraints one need to fix the complete 
particle content at high energies in order to ensure the cancelation of ultraviolet 
divergencies in arbitrary scattering process. 

The relations © were derived for effective low-energy parameters accounting 
for radiation corections. Nevertheless, they also hold at tree-level in a wide class 
of known models containing the Abelian Z' . In this case, it is possible to derive 
them by imposing the requirement that the SM Lagrangian (including Yukawa 
term) has to be invariant with respect to the extra V(V) group associated with 
the Z' PJ]. 

A lot widely discussed models are derived from the Eg group (the so called 
LR, X"V" models). The tree- level Z' couplings to the SM fermions in the models 
are shown in Table [TJ 
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The Eg-symmetry breaking scheme 

E 6 -> SO(10) x U(1V, SO(10) -> SU(3) C x SU(2) L x SU(2) fi x U(1) B _ £ . 

leads to the so called left-right (LR) model. Another scheme, 

E 6 -> SO(10) x U(l) v > -> SU(5) x U(l) x x U(l)^, 

predicts the Abelian Z' which is a linear combination of the neutral vector 
bosons ip and x> 

Z' = \ cos (3 + ip sin (3 

with the mixing angle /3. If we suppose only one Z' boson at low energies, the 
ip boson should be much heavier than the \ held. In this case the field ijj is 
decoupled and j3 0. As it is seen, both the LR and the X'i 3 models (with 
f3 = to avoid two Z' bosons with the same scale of masses) satisfy the relations 
(JTJ) except for neutrinos. This fact is a consequence of the zero neutrino mass 
assumed already. Neutrinos are not detected in the discussed experiments, so 
the question about Z 1 interactions to neutrinos is inessential. 

The relations ([7]) are valid not only for the E§ based models in Table [TJ For 
example, the relations also cover the Sequential SM (SSM) mentioned in reports 
of LEP Collaborations. Thus, they describe correlations between Z' couplings 
for a wide set of models beyond the SM. That is the reason to call the relations 
model-independent ones. 

LEP collaborations have applied model dependent search for Z 1 and obtained 
the low bounds on the mass my > 400 — 800 GeV dependently on a specific 
model [TJ [2J [31 HI [5] . In our analysis, the relations ([TJ give a possibility to reduce 
the number of fitted parameters, to determine kinematics of the processes, and 
to introduce observables which uniquely pick out the Z' signals. Therefore we 
are able to distinguish the particle instead of constraining its mass. 

3 Z' search in e + e — > [i + [i~,t + t~ processes 
3.1 The differential cross section 

Let us consider the processes e + e~ — » (I = h,t) with the non-polarized 

initial and final state fermions. In order to introduce the observable which 
selects the signal of the Abelian Z' boson we need to compute the differential 
cross-sections of the processes up to the one-loop level. Two classes of Feynman 
diagrams are taken into account. The first one includes the pure SM graphs. The 
set of SM diagrams give the SM prediction for the process which is a background 
for observation of possible deviations due to Z 1 boson. Obviously, the SM has 
to be estimated as accurate as possible. So, the full set of radiative corrections 
must be taken into account. They are the mass operators, the vertex corrections, 
the boxes, and the effects of initial and final state radiation of soft photons. The 
kinematic region allowed by the detectors is also important. Fine cancelations 
of ultraviolet and infrared divergencies occur due to renormalizability of the SM 
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giving the finite result. The resulting SM value is published by the Collaboration 
and have been checked. 

The second class of diagrams includes heavy Z' boson as a virtual state. 
Such graphs lead to small corrections of order s/m 2 z , to cross section. Since 
the effective low energy Lagrangian is used to describe Z' interactions to SM 
particles and the particle content of the underlying theory remains hidden, one 
has to consider the Z' contribution in the decoupling limit. Namely, we assume 
that Z' is not excited inside loops. The tree-level diagram e + e~ — > Z' — > 
defines a leading contribution to the cross-section. It is enough to take 
into account this diagram to estimate possible Z' signals. The cross-section 
includes the interference of the Z'-exchange amplitude with the SM amplitudes. 
Radiative corrections were incorporated with the Z'-exchange diagram in the 
improved Born approximation. 

In actual calculations and experimental data treating, the SM values of cross 
sections coincide with the results of the LEP Collaborations and the deviations 
due to the Z' boson have been computed in the improved Born approximation 
at one-loop level. The same approach is used for the Bhabha process which 
will be analyzed in next section. This is sufficient to analyze the present day 
experimental data. 

In the lower order in m^, 2 the Z' contributions to the differential cross-section 
of the process e + e~ — > l + l~ are expressed in terms of four-fermion contact cou- 
plings, only. If one takes into consideration the higher-order corrections in m^ 2 , 
it becomes possible to estimate separately the Z'-induced contact couplings and 
the Z' mass [21]. In the present analysis we keep the terms of order 0{m z f) to 
fit both of these parameters. 

Expanding the differential cross-section in the inverse Z' mass and neglecting 
the terms of order 0(m z f), we have 

-^r 1 +EE l A \M *) + 4-(*> *)c] **** 

Z i=l 3=1 

7 i j k 

+ zZ/ZzZzZ C ^kn( S ^ Z )^ X J X kXn, (9) 
z—l j — 1 k—1 n—1 

where the dimensionless quantities 

2 

C=— (X 1 ,X2,X 3 ,X4,,X5,X 6 ,X 7 ) = (a,V e ,Va,V T ,V d ,V s ,V b ) (10) 

are introduced. Since the axial-vector couplings of the Abelian Z' boson are 
universal, we use the shorthand notation a = a e . In what follows the index 
I = fi,r denotes the final-state lepton. 

The coefficients A, B, C are determined by the SM couplings and masses. 
Each factor may include the tree-level contribution, the one-loop correction 
and the term describing the soft-photon emission. The factors A describe the 
leading-order contribution, whereas others correspond to the higher order cor- 
rections in m z , 2 . 
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3.2 The observable 

To take into consideration the correlations ([3]) we introduce the observable (Ji(z) 
defined as the difference of cross sections integrated in some ranges of the scat- 
tering angle 8 [T31 [TS] : 

ai(z) = [ dai dcosO - [ d<Jl rfcosfl, (11) 



dcos8 dcos 

where z stands for the cosine of the boundary angle. The idea of introducing the 
z-dependent observable (TTT1) is to choose the value of the kinematic parameter 
z in such a way that to pick up the characteristic features of the Abelian Z' 
signals. 

The deviation of the observable from its SM value can be derived by the 
angular integration of the differential cross-section and has the form: 



Act,(z) = a l (z)-af u (z)=J2 1 £[ A U s ^)+ & iM^ 

i=l j=l 
7 i j k 

+J2J2J2J2^kn(^ z ) x -^j x kx n . (12) 

i— 1 j=l k=l n— 1 

Then let us introduce the quantity Act (z) = ct (z) — <jsm (z) which owing to 
the relations ([7]) can be written in the form 

Aaj(z) = ^^[Fi(z,s)Y*+2F((z, S )T 3f Y L , f Y L , e 

+ 2Fi(z,s)T 3f Y L jY$ + Fi(z,s)Y L , e Y^ . (13) 

The factor functions F/(z,s) depend on the fermion type through the \Qf\, 
only. In Fig. [T]they are shown as the functions of z for ^/s = 500 GeV. The 
leading contributions to F/(z,s), 

Fi(z,s) = -5lQ/l(i---- 2 -y)+°(^), 

F{(z,s) = t[i-^-\Q f \(sz + ^)]+o(^y 

F*{z,s) = -\ ^-z 2 ) + \ (3^ + ^ 3 ) (4|Q/|-1) + o(- 

Fi(z,s) = l\Q f \(l-3z-z*-z 3 )+o(^fy (14) 

are given by the Z' exchange diagram e~e + — > Z' — > //, since the Z—Z' mixing 
contribution to the Z exchange diagram is suppressed by the factor m z /s. 
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Table 1: The Z' couplings to the SM fermions in the most discussed E 6 -based 



models. 
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Figure 1: The leptonic functions Fq (the solid curve), F[ (the long-dashed 
curve), F% (the dashed curve), and (the dotted curve) at y/s = 500 GcV. 
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From Eqs. (fl"4")l one can see that the leading contributions to the leptonic 
factors F[, F%, are found to be proportional to the same polynomial in z. 
This is the characteristic feature of the leptonic functions F\ originating due 
to the kinematic properties of fermionic currents and the specific values of the 
SM leptonic charges. Therefore, it is possible to choose the value of z = z* 
which switches off three leptonic factors F[, F%, F^ simultaneously. Moreover, 
the quark function F^ in the lower order is proportional to the leptonic one and 
therefore is switched off, too. As is seen from Fig. [TJ the appropriate value of 
z* is about ~ 0.3. By choosing this value of z* one can simplify Eq. (fT3|). It 
is also follows from Eq. (flU)) that neglecting the factors F[, F^, F^ one obtains 
the sign definite quantity Aai(z*) ~ Y? ~ a 2 . 

There is the interval of boundary angle values at which the factors A\ 1: B l n , 
and C l nll at the sign-definite parameters a 2 , a 2 £, and a 4 contribute more than 
95% of the observable value. It gives a possibility to construct the sign-definite 
observable Acr;(z*) < by specifying the proper value of z* . 

In general, one could choose the boundary angle z* in different schemes. If 
just a few number of tree-level four-fermion contact couplings are considered, 
one can specify z* in order to cancel the factor at the vector-vector coupling. 
However, if one-loop corrections are taken into account there is a large amount of 
additional contact couplings. So, we have to define some quantitative criterion 
F(z) to estimate the contributions from sign-definite factors at a given value of 
the boundary angle z. Maximizing the criterion, one could derive the value z* 
corresponding to the sign-definite observable Acr;(z*). Since the observable is 
linear in the coefficients A, B, and C, we introduce the following criterion, 

(15) 
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where the positive 'weights' ~ ( and ojc ~ e take into account the order of 
each term in the inverse Z' mass. 

The numeric values of the 'weights' ujb and ujq can be taken from the present 
day bounds on the contact couplings pQ. As the computation shows, the value 
of z* with the accuracy 10~ 3 depends on the order of the 'weight' magnitudes, 
only. So, in what follows we take wb~4x 10 -3 and wc ~4x 10~ 5 . 

The function z*(s) is the decreasing function of the center-of-mass energy. 
It is tabulated for the LEP2 energies in Table [2] The corresponding values of 
the maximized function F are within the interval 0.95 < F < 0.96. 

Since A[ 1 (s, z*) < 0, B l n (s,z*) < and C l nn (s,z*) < 0, the observable 

Aa t (z*) = [A l n (s, z*) + (& n (s, z*)] a 2 + C l lnl (s, z*)a 4 (16) 

is negative with the accuracy 4-5%. Since this property follows from the rela- 
tions ([5]) for the Abelian Z' boson, the observable Aai(z*) selects the model- 
independent signal of this particle in the processes e + e~ — > . It allows to 
use the data on scattering into \i\x and rr pairs in order to estimate the Abelian 
Z' coupling to the axial- vector lepton currents. 
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Although the observable can be computed from the differential cross-sections 
directly, it is also possible to recalculate it from the total cross-sections and the 
forward-backward asymmetries. The recalculation procedure has the proper 
theoretical accuracy. Nevertheless, it allows to reduce the experimental errors 
on the observable, since the published data on the total cross-sections and the 
forward-backward asymmetries are more precise than the data on the differential 
cross-sections. 

The recalculation is based on the fact that the differential cross-section can 
be approximated with a good accuracy by the two-parametric polynomial in the 
cosine of the scattering angle z: 

do* (a) daf M (s) M 2 . 

— = i + C 1 + + zm + 8i(z), 17 
dz dz 

where 5i(z) measures the difference between the exact and the approximated 
cross-sections. The approximated cross-section reproduces the exact one in the 
limit of the massless initial- and final-state leptons and if one neglects the con- 
tributions of the box diagrams. Detailed analysis of this point is given in [15] 
where it was showed that the the theoretical error is one order less than the 
corresponding statistical uncertainty for the observable. Thus, the proposed 
approximation is quite good and can be successfully used to obtain more accu- 
rate experimental values of the observable. 



3.3 Data fit 

To search for the model-independent signals of the Abelian Z'-boson we will 
analyze the introduced observable Aai(z*) (1 16|) on the base of the LEP2 data 
set. In the lower order in m^ 2 it depends on one flavor- independent parameter 
a 2 , 

Aa^{z*) = i' u ( S , z*)d 2 + & llu (a, z*)d\ (18) 

which can be fitted from the experimental values of Aa fl (z*) and Aa T (z*). As 
we axplained above, the sign of the fitted parameter (a 2 > 0) is the characteristic 
feature of the Abelian Z' signal. 

In what follows we will apply the usual fit method based on the likelihood 
function. The central value of d 2 is obtained by the minimization of the % 2 - 
function: 



x 2 (d 2 ) = y [^n(*-)-A^V)r (19) 



where the sum runs over the experimental points entering a data set chosen. 
The la CL interval (&i, b-i) for the fitted parameter is derived by means of the 
likelihood function C(d 2 ) oc exp[— % 2 (a 2 )/2]. It is determined by the equations: 

/ 2 C(e')de' = 0.68, C(h) = £(b 2 ). (20) 
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To relate our results with those of Refs. [I] we introduce the contact inter- 
action scale 

A 2 = 4m|a- 2 . (21) 

This normalization of contact couplings is admitted in Refs. !lj. We use again 
the likelihood method to determine a one-sided lower limit on the scale A at 
the 95% CL. It is derived by the integration of the likelihood function over the 
physically allowed region a 2 > 0. The strict definition is 



y>e* poo 

A = 2m z (e*)- 1 / 2 , / C(e')de' = 0.95 / C(e')de'. 

Jo Jo 



(22) 



We also introduce the probability of the Abelian Z' signal as the integral of 
the likelihood function over the positive values of a 2 : 



P 



L(e')de'. 



(23) 



Actually, the fitted value of the contact coupling a 2 originates mainly from 
the leading-order term in the inverse Z' mass in Eq. (jTHJ) - The analysis of 
the higher-order terms allows to estimate the constraints on the Z' mass alone. 
Substituting a 2 in the observable (fTS)) by its fitted central value, one obtains 
the expression 



Aa l { Z *) = ii^ (*,«*)+ C^il 4ttcd + ^ m (^*)«fitted 



(24) 



which depends on the parameter £ — m 2 z /m 2 z , . Then, the central value of this 
parameter and the corresponding la CL interval are derived in the same way 
as those for a 2 . 

To fit the parameters a 2 and £ we start with the LEP2 data on the total 
cross-sections and the forward-backward asymmetries pQ. The corresponding 
values of the observable Aai(z*) with their uncertainties 5ai(z*) are calculated 
from the data by means of the following relations: 



6<ri(z*) 2 



t( 3 +^ 2 ) 
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{l-z* 2 )a^ SM AAf B , 



A r(l-z* 2 )-^(3 + z* 2 ) 
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We perform the fits assuming several data sets, including the fj.fi, tt, and 
the complete fifi and tt data, respectively. The results are presented in Table El 
As is seen, the more precise fifi data demonstrate the signal of about la level. 
It corresponds to the Abelian Z'-boson with the mass of order 1.2-1.5 TeV if 
one assumes the value of a = £r 2 /47r to be in the interval 0.01-0.02. No signal is 
found by the analysis of the tt cross-sections. The combined fit of the fifi and 
tt data leads to the signal below the la CL. 
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Table 2: The boundary angle z* and the coefficients in the observable Aai(z*) 
for the scattering into /i and r pairs at the one-loop level. 



v 7 ?, GeV 


2 


* 


An 


xlO 2 






T + T~ 




r + T - 


130 


0.450 


0.460 


-7.29 


-6.87 


136 


0.439 


0.442 


-7.09 


-6.88 


161 


0.400 


0.400 


-6.43 


-6.25 


172 


0.390 


0.391 


-6.19 


-6.01 


183 


0.383 


0.385 


-5.99 


-5.71 


189 


0.380 


0.380 


-5.86 


-5.68 


192 


0.380 


0.380 


-5.79 


-5.62 


196 


0.380 


0.379 


-5.71 


-5.54 


200 


0.378 


0.378 


-5.64 


-5.47 


202 


0.376 


0.377 


-5.60 


-5.43 


205 


0.374 


0.374 


-5.55 


-5.48 


207 


0.372 


0.372 


-5.52 


-5.44 



Table 3: The contact coupling a 2 with the 68% CL uncertainty, the 95% CL 
lower limit on the scale A, the probability of the Z' signal, P, and the value of 
C = m 2 z /m z , as a result of the fit of the observable recalculated from the total 
cross-sections and the forward-backward asymmetries. 



Data set 


a A , xlQ- b 


A, TcV 


P 


C, xlO" 2 




^•00-4.86 


16.4 


0.77 


0.9 ±27.8 


TT 


-2 66+ 6 43 

z - OD -6.39 


17.4 


0.34 


-0.1 ±50.1 


and tt 


-i oo+3.89 
- L - 00 -3.87 


19.7 


0.63 


1.7 ±60.9 
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Being governed by the next-to-leading contributions in m z J , the fitted values 
of £ are characterized by significant errors. The data set gives the central 
value which corresponds to mz' ~ 1.1 TeV. 

We also perform a separate fit of the parameters based on the direct cal- 
culation of the observable from the differential cross-sections. The experimen- 
tal uncertainties of the data on the differential cross-sections are of one order 
larger than the corresponding errors of the total cross-sections and the forward- 
backward asymmetries. These data also provide the larger values of the contact 
coupling a 2 . As for the more precise /x/x data, three of the LEP2 Collaborations 
demonstrate positive values of a 2 . The combined a 2 is also positive and remains 
practically unchanged by the incorporation of the rr data. 



4 Search for Z' in e + e — > e + e process 
4.1 The differential cross-section 

In our analysis of the Bhabha process, as the SM values of the cross-sections 
we use the quantities calculated by the LEP2 collaborations [3J [3J 2J [SJ [331 H3] • 
They account for either the one-loop radiative corrections or initial and final 
state radiation effects (together with the event selection rules, which are specific 
for each experiment). As it is reported by the DELPHI Collaboration, there is 
a theoretical error of the SM values of about 2%. In our analysis this error is 
added to the statistical and systematic ones for all the Collaborations. As it was 
checked, the fit results are practically insensitive to accounting for this error. 

The deviation from the SM is computed in the improved Born approxima- 
tion. This approximation is sufficient for our analysis leading to the systematic 
error of the fit results less than 5-10 per cents. 

The deviation from the SM of the differential cross-section for the process 
e + e~ — > £ + £~ can be expressed through various quadratic combinations of cou- 
plings a — a e , v e , v^, v T . For the Bhabha process it reads 

= ^ + ^ + <*> 

where the factors are known functions of the center-of-mass energy and the 
cosine of the electron scattering angle z plotted in Fig. [5] The deviation of the 
cross-section for e + e~ — > /i + fjT (t + t~) processes has a similar form 

az az m z , m z , m z , m z , 

Eqs. ([2g ]) -((2T |) are our definition of the Z' signal. 

Note again that the cross-sections in Eqs. (|26[) -(|27 |) account for the relations 
10 through the functions fi(z), f^{z), fi(z), since the coupling (the mixing 
angle 9q) is substituted by the axial coupling constant a. Usually, when a four- 
fcrmion effective Lagrangian is applied to describe physics beyond the SM [53] , 
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Figure 2: The factors at the Z' couplings in the differential cross-section of the 
Bhabha process. 




Figure 3: Factors F a (y/s,z) (solid) and F av (* s fs,z) (dashed) in the normalized 
deviation of the differential cross-section da / dz for yfs = 200 and 500 GeV. 



this dependence on the scalar field coupling is neglected at all. However, in our 
analysis when we are interested in searching for signals of the Z'-boson on the 
base of the effective low-energy Lagrangian ©-(H}, these contributions to the 
cross-section are essential. 



4.2 One-parameter fit 

The factor /f e (z) is positive monotonic function of z (see Fig. [2] for the center- 
of-mass energies y/s = 200 GeV. The same behavior is observed for higher ener- 
gies). Such a property allows one to choose fi e (z) as a normalization factor for 
the differential cross section. Then the normalized deviation of the differential 
cross-section reads [16] 

T- = A 7ee(, A ^T=v 2 + F a (v/i, z)a 2 + F OT (v/i, z)av + ..., (28) 
dz 4nj2 (z) dz 
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and the normalized factors are shown in Fig [3] for energies of LEP and ILC 
experiments. Now these factors are finite at z — > 1. Each of them in a special 
way influences the differential cross-section. 

1. The factor at v 2 is just the unity. Hence, the four-fermion contact coupling 
between vector currents, v 2 , determines the level of the deviation from the 
SM value. 

2. The factor at a 2 depends on the scattering angle in a non-trivial way. It 
allows to recognize the Abelian Z' boson, if the experimental accuracy is 
sufficient. 

3. The factor at av results in small corrections. 

Thus, effectively, the obtained normalized differential cross-section is a two- 
parametric function. In the next sections we introduce the observables to fit 
separately each of these parameters. 

4.3 Observables to pick out v 2 

The normalized deviation of the differential cross-section (|2"5)) is (effectively) 
the function of two parameters, a 2 and v 2 . We are going to introduce the 
integrated observables which determine separately the four-fermion couplings 
a 2 and v 2 [TB] . 

Let us first proceed with the observable for v 2 . After normalization the 
factor at the vector-vector four-fermion coupling becomes the unity. Whereas 
the factor at a 2 is a sign-varying function of the cosine of the scattering angle. 
As it follows from Fig. [3J for the center-of-mass energy 200 GeV it is small over 
the backward scattering angles. So, to measure the value of v 2 the normalized 
deviation of the differential cross-section has to be integrated over the backward 
angles. For the center-of-mass energy 500 GeV the factor at a 2 is already a non- 
vanishing quantity for the backward scattering angles. The curves corresponding 
to intermediate energies are distributed in between two these curves. Since they 
are sign-varying ones at each energy point some interval of z can be chosen to 
make the integral to be zero. Thus, to measure the Z' coupling to the electron 
vector current v 2 we introduce the integrated cross-section (|28|) 



where at each energy the most effective interval [zq, zq + Az] is determined by 
the following requirements: 

1. The relative contribution of the coupling v 2 is maximal. Equivalently, the 
contribution of the factor at a 2 is suppressed. 

2. The length Az of the interval is maximal. This condition ensures that the 
largest number of bins is taken into consideration. 




(29) 



1G 





-0.6 -0.4 -0.2 z n 



-0.4 -0.2 z n 



Figure 4: Relative contribution of the factor at v 2 to the observable cry as the 
function of the left boundary of the angle interval, zq, and the interval length, 
Az, at the center-of-mass energy 200 and 500 GeV. The shaded areas correspond 
to the contributions > 95% (dark), from 90% to 95% (midtone), and < 90% 
(light). 



The relative contribution of the factor at v 2 is defined as 

Az 

Ky = 



zo+Az 



F a dz 



z +Az 



F dz 



(30) 



and shown in Fig. |4] as the function of the left boundary of the angle inter- 
val, zq, and the interval length, Az. In each plot the dark area corresponds 
to the observables which values are determined by the vector-vector coupling 
v 2 with the accuracy > 95%. The area reflects the correlation of the width 
of the integration interval Az with the choice of the initial zq following from 
the mentioned requirements. Within this area we choose the observable which 
includes the largest number of bins (largest Az). The corresponding values of 
Zq and Az are marked by the white dot on the plots in Fig. 0J As the carried 
out analysis showed, the point zq is shifted to the right with increase in energy 
whereas Az remains approximately the same. 

From the plots it follows that the most efficient intervals are 



-0.6 < z < 0.2, y/s = 200 GeV, 
-0.3 < z < 0.7, \/s = 500 GeV. 



(31) 



Therefore the observable (|2"9l allows to measure the Z' coupling to the electron 
vector current v 2 with the efficiency > 95%. 

Fitting the LEP2 final data with the one-parameter observable, we find the 
values of the Z' coupling to the electron vector current together with their la 
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uncertainties: 



ALEPH 


«£ = 


-0.11 ±6.53 x 10~ 4 


DELPHI 


Ve = 


1.60 ± 1.46 x 10~ 4 


L3 


vl = 


5.42 ±3.72 x 10~ 4 


OPAL 


vl = 


2.42 ± 1.27 x 10~ 4 


Combined 


vl = 


2.24 ±0.92 x 10~ 4 . 



As one can see, the most precise data of DELPHI and OPAL collaborations 
are resulted in the Abelian Z' hints at one and two standard deviation level, 
correspondingly. The combined value shows the 2a hint which corresponds to 
0.006 < \v e \ < 0.020. 



4.4 Observables to pick out a 2 

In order to pick the axial- vector coupling a 2 one needs to eliminate the domi- 
nant contribution coming from v 2 . Since the factor at v 2 in the da/dz equals 
unity, this can be done by summing up equal number of bins with positive and 
negative weights. In particular, the forward-backward normalized deviation of 
the differential cross-section appears to be sensitive mainly to a 2 , 

<7fb = / ^ Z ~dz \ dz —F ayFB d 2 + F aVyFB dv. (32) 

The value z max is determined by the number of bins included and, in fact, 
depends on the data set considered. The LEP2 experiment accepted e+e~ 
events with \z\ < 0.72. In what follows we take the angular cut z max = 0.7 for 
defmitcness. 

The efficiency of the observable is determined as: 

I^fbI (33) 



|-Fa,FB| + l-Fa?, FB| 



It can be estimated as k = 0.9028 for the center-of-mass energy 200 GeV and 
k = 0.9587 for 500 GeV. Thus, the observable 

CT FB = 0.224a 2 - 0.024aw, V« = 200 GeV, 

ctfb = 0.472a 2 - 0.020aw, sfs = 500 GeV (34) 

is mainly sensitive to the Z' coupling to the axial- vector current a 2 . 

Consider a usual situation when experiment is not able to recognize the an- 
gular dependence of the differential cross-section deviation from its SM value 
with the proper accuracy because of loss of statistics. Nevertheless, a unique 
signal of the Abelian Z' boson can be determined. For this purpose the observ- 
ables fz° +Az {da I 'dz)dz and <tfb must be measured. Actually, they are derived 
from the normalized deviation of the differential cross-section. If the deviation 
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is inspired by the Abelian Z' boson both the observables are to be positive 
quantities simultaneously. This feature serves as the distinguishable signal of 
the Abelian Z' virtual state in the Bhabha process for the LEP2 energies as well 
as for the energies of future electron-positron collider ILC (> 500 GeV). The 
observables fix the unknown low energy vector and axial- vector Z' couplings to 
the electron current. Their values have to be correlated with the bounds on a 2 
and v 2 derived by means of independent fits for other scattering processes. 

We estimated the observable related to the value of a 2 . Since in the 
Bhabha process the effects of the axial-vector coupling are suppressed with 
respect to those of the vector coupling, we expect much larger experimental 
uncertainties for a 2 . Indeed, the LEP2 data lead to the huge errors for a 2 of 
order 10 -3 — 10~ 4 . The mean values are negative numbers which are too large 
to be interpreted as a manifestation of some heavy virtual state beyond the 
energy scale of the SM. 

Thus, the LEP2 data constrain the value of v 2 at the 2a CL which could cor- 
respond to the Abelian Z' boson with the mass of the order 1 TeV. In contrast, 
the value of a 2 is a large negative number with a significant experimental un- 
certainty. This can not be interpreted as a manifestation of some heavy virtual 
state beyond the energy scale of the SM. 



4.5 Many-parameter fits 



To account for all the accessible data of LEP experiments, we address to many 
parameter fits [T7]. As the basic observable to fit the LEP2 experiment data on 
the Bhabha process we propose the differential cross-section 



da 



Bhabha 



da 



Bhabha.SAf 



dz 



dz 



(35) 



where i runs over the bins at various center-of-mass energies ^fs. The final 
differential cross-sections measured by the ALEPH [22] (130-183 GeV), DELPHI 
(189-207 GeV), L3 [23] (183-189 GeV), and OPAL [2][3]|4] (130-207 GeV) 
collaborations are taken into consideration (299 bins). 



As the observables for e + e 



A t+ A t ~ 



processes, we consider the total 



cross-section and the forward-backward asymmetry 



(Jrp 



i + r ,sm 



.4 



FB 



A 



rr,sM 

FB 



(36) 



where i runs over 12 center-of-mass energies y/s from 130 to 207 GeV. We 
consider the combined LEP2 data pQ for these observables (24 data entries 
for each process). These data are more precise as the corresponding differen- 
tial cross-sections. Our analysis is based on the fact that the kinematics of 
s-channel processes is rather simple and the differential cross-section is effec- 
tively a two-parametric function of the scattering angle. The total cross-section 
and the forward-backward asymmetry incorporate complete information about 
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the kinematics of the process and therefore are an adequate alternative for the 
differential cross-sections. 

The data are analysed by means of the x 2 fit [T7] . Denoting the observables 
(f55|) — by (7j, one can construct the x 2 -function, 

X 2 (a,v e ,v^,v T ) = ^2 

i 

where cr ox and 8a are the experimental values and the uncertainties of the ob- 
servables, and cr th are their theoretical expressions presented in Eqs. (|26p - (j2"?| . 
The sum in Eq. (j37]) refers to either the data for one specific process or the com- 
bined data for several processes. By minimizing the x 2 -function, the maximal- 
likelihood estimate for the Z 1 couplings can be derived. The x 2 -function is also 
used to plot the confidence area in the space of parameters a, v e , v^, and v T . 
Note that in this way of experimental data treating all the possible correla- 
tions are neglected. We believe that at the present stage of investigation this is 
reasonable, because the Collaborations have never reported on this possibility. 

For all the considered processes, the theoretic predictions af 1 are linear com- 
binations of products of two Z' couplings 

7 

af 1 = y~]CjjAj, Aj = {a?,Ve,a%,%v^,%v r ,aVii,a,v r }, (38) 
i=i 

where CV,- are known numbers. 

In the Bhabha process, the Z' effects are determined by three linear-independent 
contributions coming from a 2 , ?) 2 , and av e and the number degrees of freedom 
(d.o.f) M — 3. As for the e + e~ —> /i + fi~ ', t 4 't~ processes, the observables de- 
pend on four linear-independent terms for each process: a 2 , v e v^, v e a, av^ for 
e + e~ — > i-i + fJ,~; and a 2 , v e v T , v e a, av T for e + e _ — » r + r _ (M = 4). Note that 
some terms in the observables for different processes are the same. Therefore, 
the number of \ 2 d.o.f. in the combined fits is less than the sum of d.o.f. for 
separate processes. Hence, the predictive power of the larger set of data is not 
drastically spoiled by the increased number of d.o.f. In fact, combining the data 
of the Bhabha and e + e~ — > (t + t~) processes together we have to treat 

five linear-independent terms. The complete data set for all the lepton processes 
is ruled by seven d.o.f. As a consequence, the combination of the data for all 
the lepton processes is possible. 

The parametric space of couplings (a, v e , w M , v T ) is four-dimensional. How- 
ever, for the Bhabha process it is reduced to the plane (a, v e ) 7 and to the 
three-dimensional volumes (a, v e , v^), (a, v e , v T ) for the e + e~ — > and 
e + e~ — > t + t~ processes, correspondingly. The predictive power of data is 
distributed not uniformly over the parameters. The parameters a and v e are 
present in all the considered processes and appear to be significantly con- 
strained. The couplings or v T enter when the processes e + e~ — > fi fi~ 
or e + e~ — > t + t~ are accounted for. So, in these processes, we also study the 
projection of the confidence area onto the plane (a,v e ). 



5a, 



(37) 
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Figure 5: The 95% CL areas in 

the (a, v e ) plane for the Bhabha, 

e + e~ — > and e + e~ — » t + t~ 
processes. 



Figure 6: The projection of the 
95% CL area onto the (a, v e ) plane 
for the combination of the Bhabha, 

e + e~ — » and e + e~ — > t + t~ 

processes. 



The origin of the parametric space, a — v e = 0, corresponds to the absence 
of the Z' signal. This is the SM value of the observables. This point could 
occur inside or outside of the confidence area at a fixed CL. When it lays out 
of the confidence area, this means the distinct signal of the Abelian Z'. Then 
the signal probability can be defined as the probability that the data agree with 
the Abelian Z' boson existence and exclude the SM value. This probability 
corresponds to the most stringent CL (the largest Xcl) at which the point 
a — v e — is excluded. If the SM value is inside the confidence area, the Z' 
boson is indistinguishable from the SM. In this case, upper bounds on the Z' 
couplings can be determined. 

The 95% CL areas in the (a, v e ) plane for the separate processes are plotted 
in Fig. [5] As it is seen, the Bhabha process constrains both the axial-vector 
and vector couplings. As for the e + e~ — > f-i + ^~ and e+e~ — > t + t~ processes, 
the axial-vector coupling is significantly constrained, only. The confidence areas 
include the SM point at the meaningful CLs, so the experiment could not pick 
out clearly the Abelian Z' signal from the SM. An important conclusion from 
these plots is that the experiment significantly constrains only the couplings 
entering sign-definite terms in the cross-sections. 

The combination of all the lepton processes is presented in Fig. [S] There is 
no visible signal beyond the SM. The couplings to the vector and axial-vector 
electron currents are constrained by the many-parameter fit as \v e \ < 0.013, 
\a\ < 0.019 at the 95% CL. If the charge corresponding to the Z' interactions 
is assumed to be of order of the electromagnetic one, then the Z 1 mass should 
be greater than 0.67 TeV. For the charge of order of the SM SU(2)l coupling 
constant mz' > 1.4 TeV. One can see that the constraint is not too severe to 
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Figure 7: The 68% CL area in the (a, v e ) plane from the backward bins of the 
Bhabha process in the LEP2 experiments (the shaded area). The hatched area 
is the 68% CL area from the LEP 1 data on the Bhabha process. 

exclude the Z' searches at the LHC. 

Let us compare the obtained results with the one-parameter fits. As one can 
see, the most precise data of DELPHI and OPAL collaborations are resulted in 
the Abelian Z' hints at one and two standard deviation level, correspondingly. 
The combined value shows the 2a hint, which corresponds to 0.006 < \v e \ < 
0.020. On the other hand, our many-parameter fit constrains the Z' coupling 
to the electron vector current as \v e \ < 0.013 with no evident signal. Why does 
the one-parameter fit of the Bhabha process show the 2a CL hint whereas there 
is no signal in the two-parameter one? Our one-parameter observable accounts 
mainly for the backward bins. This is in accordance with the kinematic features 
of the process: the backward bins depend mainly on the vector coupling v\, 
whereas the contributions of other couplings are kinematically suppressed (see 
Fig. [2]). Therefore, the difference of the results can be inspired by the data sets 
used. To clarify this point, we perform the many-parameter fit with the 113 
backward bins (z < 0), only. The \ 2 minimum, Xmm = 93.0, is found in the 
non-zero point \a\ = 0.0005, v e = 0.015. This value of the Z' coupling v e is in 
an excellent agreement with the mean value obtained in the one-parameter fit. 
The 68% confidence area in the (a,v e ) plane is plotted in Fig. [7] There is a 
visible hint of the Abelian Z' boson. The zero point a = v e = (the absence 
of the Z 1 boson) corresponds to x 2 = 97.7. It is covered by the confidence area 
with 1.3a CL. Thus, the backward bins show the 1.3cr hint of the Abelian Z' 
boson in the many-parameter fit. So, the many-parameter fit is less precise than 
the analysis of the one-parameter observables. 

At LEP1 experiments [26] the Z-boson couplings to the vector and axial- 
vector lepton currents (gv, 9a) were precisely measured. The Bhabha process 
shows the la deviation from the SM values for Higgs boson masses run > 114 
GeV (see Fig. 7.3 of Ref. [IS]). This deviation could be considered as the effect 
of the Z-Z' mixing. It is interesting to estimate the bounds on the Z' couplings 
following from these experiments. 

Due to relations J5]), the Z-Z' mixing angle is completely determined by the 
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Table 4: The summary of the fits of the LEP data for the dimensionless contact 
couplings (j6]). 



Data v'i, xlCT 4 


a 2 , x!CT b 


LEP1 


e~e+, 68% CL 


1.25 ± 1.25 


LEP2, one-parameter fits 


e~e+, 68% CL 2.24 ±0.92 




fj.fi, 68% CL 
fifi.TT, 68% CL 


RR+4.89 
O.DD_4 gg 

1 oq+3.89 
i - OJ -3.87 


LEP2, many-parameter fits 


e~e + , fifi,TT, 95% CL < 1.69 


< 36.1 


e-e+ backward, 68% CL 2.25^11% 


< 94.9 



axial- vector coupling a. So, the deviations of gy, g A from their SM values are 
governed by the couplings a and v e , 

9V ~ gv M = -49.06a« e , g A - gf 1 = 49.06a 2 . (39) 

Let us assume that the total deviation of theory from experiments follows due 
to the Z-Z' mixing. This gives an upper bound on the Z' couplings. In this 
way one can estimate whether the Z' boson is excluded by the experiments or 
not. 

The la CL area for the Bhabha process from Ref. [35] is converted into the 
(a, Ve) plane in Fig. [7J The SM values of the couplings correspond to the top 
quark mass m t — 178 GeV and the Higgs scalar mass run = 114 GeV. As it is 
seen, the LEP1 data on the Bhabha process is compatible with the Abelian Z' 
existence at the la CL. The axial-vector coupling is constrained as \a\ < 0.005. 
This bound corresponds to a 2 < 2.5 x 10~ 5 , which agrees with the one-parameter 
fits of the LEP2 data for e + e~ -> n + fi~ , t + t~ processes (a 2 = 1.3±3.89 x 10~ 5 
at 68% CL). On the other hand, the vector coupling constant v e is practically 
unconstrained by the LEP1 experiments. 

For the convenience, in Table [4] we collect the summary of the fits of the 
LEP data in terms of dimensionless contact couplings ([5]). From the analysis 
carried out we come to conclusion that, in principle, the LEP experiments were 
able to detect the Z'-boson signals if the statistics had been sufficient. 

5 Model independent results and search for Z' 
at the LHC 

In this section we discuss all the assumptions giving a possibility to pick out the 
Z' signal and determine its characteristics in a model independent way. We also 
note the role of the present results for the LHC and future ILC experiments. 

As it was already noted, in searching for this particle at the LEP and Teva- 
tron a model dependent analysis was applied. As the main motivation for this 
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approach it was the different number of chiral fermions involved in different 
models (see, for example, Ref. [7]). In this way the low bounds on mz* have 
been estimated and the smallness of the Z-Z' mixing was also observed. 

On the contrary, in our approach the relations between the parameters of 
the effective low energy Lagrangians have been accounted for that gave a possi- 
bility to determine not only the bounds but also the mass and other parameters 
of the Z'. 

To be precise, let us note all the assumptions used in our investigations. We 
analyzed the four-fermion scattering amplitudes of order ~ m^, generated by 
the Z' virtual states. The vertices linear in Z' were included into the effective 
low-energy Lagrangian. We also impose a number of natural conditions. The 
interactions of a renormalizable type are dominant at low energies ~ niw- The 
non-renormalizable interactions generated at high energies due to radiation cor- 
rections are suppressed by the inverse heavy mass ~ 1 /mz> and neglected. We 
also assumed that the SU(2)l x U(1)y gauge group of the SM is a subgroup 
of the GUT group. As a consequence, all the structure constants connecting 
two SM gauge bosons with Z' have to be zero. Hence, the interactions of gauge 
fields of the types Z'W + W~ , Z' ' ZZ , and other are absent at a tree level. Our 
effective Lagrangian is also consistent with the absence of the tree-level flavor- 
changing neutral currents (FCNCs) in the fermion sector. The renormalizable 
interactions of fermions and scalars are described by the Yukawa Lagrangian 
that accounts for different possibilities of the Yukawa sector without the tree- 
level FCNCs. These assumptions are quite general and satisfied in a wide class 
of Eq inspired models. 

Within these constraints for the low energy effective Lagrangian the relations 
0,® have been derived. Correspondingly, the model independent estimates 
of the mass mz> and other parameters are regulated by the noted requirements. 
Therefore, the extended underlying model has also to accept them. 

In this regard, let us discuss the role of the obtained estimates for the LHC. 
As it is well known (see, for example, |5]), there are many tools at the LHC 
for Z' identification. But many of them are only applicable if this particle is 
relatively light. Our results are in favor to this case. 

Next important point is the determination of Z' couplings to the various 
SM fermions. As we have shown, the axial-vector couplings of the Z' to the 
SM fermions are universal and proportional to its coupling to the Higgs field. 
Hence we have obtained an estimate of the a = a/ couplings for both leptons 
and quarks. This is an essential input because experimental analysis for the 
LHC have mainly concentrated on being able to distinguish models and not on 
actual couplings^] The vector coupling v e was also estimated that, in particular, 
may help to distinguish the decay of the Z' resonance state to e + e~ pairs. Since 
the couplings a e and v e were estimated there is a possibility to distinguish this 
process from the decay of the KK system. In the literature on searching for 
the Z' it is also mentioned [5J [S7J that the determination of the Z' couplings 
to fermions could be fulfilled channel by channel, a e ,v e ,v e j,,a e ,bi ■■■■ In that 

1 Discussion of the determination of couplings can be found in Ref. I29| . 
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considerations the relations between the parameters have not been taken into 
account. But this is very essential for treating of experimental data and intro- 
ducing the relevant observables to measure. Our consideration could be useful 
in this problem. 

Other parameter is the Z—Z 1 mixing which is responsible for the different 
decay processes and the effective interaction vertices generated at the LHC 
[3 [8]. It is also determined by the axial- vector coupling (see Eq. ((8])) and 
estimated in a model independent way. Remind that in our analysis the mixing 
was systematically accounted for. Its value is of the same order of magnitude as 
the parameters that were fitted in experiments. It worth to note also that the 
existence of other heavy particles with masses mx > rnz> does not influence 
the RG relations. 

An important role of the model independent results for searching for Z' at 
the Tevatron, LHC and ILC consists, in particular, in a possibility to determine 
the particle as a virtual state due to a large amount of relevant events. We 
mentioned already that, in principle, LEP2 experiments were able to determine 
it if the statistics was sufficiently large. Experiments at the ILC will increase in 
many times the data set of interest. In fact, the observables, introduced in sects. 
6 and 7 for picking out uniquely and v\ couplings in the leptonic scattering 
process, are also effective at energies \fs > 500 GeV and could be applied in 
future experiments at ILC. 

Other model independent methods of searching for the Z' as a resonance 
state are proposed in the literature (see Refs. [27J [28l [29]). We do not discuss 
them here because they take into consideration no relations between the param- 
eters. As we mentioned already, the main goal of the present paper is to adduce 
model independent information about the Z' followed from experiments at low 
energies. Different aspects of Z' physics at the LHC are out of the scope of it. 

6 Discussion 

In this section we collect in a convenient form all the results obtained and make 
a comparison with other investigations on searching for Z' at low energies. In 
fact, this is a large area to discuss. References to numerous results obtained 
in either model dependent or model independent approaches can be found in 
Refs. [Tl [5] - Further subdivision can be done into the considerations accounting 
for any type correlations between the parameters of the low energy effective 
interactions and that of assuming complete independence of them. Because of 
a large amount of fitting parameters the latter are less predictable. 

Now, for a convenience of readers we present the results of fits of the Z' 
parameters in terms of the popular notations [6j [7] . The Lagrangian reads 

Czff = lzJrM M + Aj)- 7 5 (a s f M + Af)]f, 
Cz>ft = l Z 'Jr[(v' f -j 5 a' f )]f, (40) 
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Table 5: The summary of the fits of the LEP data for the maximum likelihood 
values of the Z 1 couplings (J40J) to the SM fermions and of the Z-Z' mixing angle 
8 . M = Y^kj denotes the unknown value of the Z 1 mass in TeV units. 



Data 


\0 O \, xHT 3 \v' e \, xlO- 1 \a' f \, xlO- 1 


Af, xlO" 3 


LEP1 


e e + 


3.17A/- 1 - 1.38M 


0.437 


LEP2, one-parameter fits 


e e + 


5.83M 






5.42M" 1 - 2.36M 


1.278 




3.27A/" 1 - 1.42M 


0.464 


LEP2, many-parameter fits 


e e+, z < 


5. 84 AT 





Table 6: The summary of the fits of the LEP data for the confidence intervals 
for the Z' couplings (l4Tfl) to the SM fermions and for the Z-Z' mixing angle 9$. 
M = pjfy denotes the unknown value of the Z' mass in TeV units. 



Data 


CL 


|0 O |, xlO" 3 \v' e \, xlO" 1 


\a' f \, xlO- 1 


A^, xlO" 3 


LEP1 


e e + 


68% 


(0;4.48)M _i 


(0; 1.95)M 


(0; 0.873) 


LEP2, one-parameter fits 


e e + 


95% 
95% 
95% 


(2.46;7.87)M 

(0; 10.39)M- 1 
(0;8.64)Af- 1 


(0;4.52)M 
(0;3.75)Af 


(0; 4.694) 
(0; 3.244) 


LEP2, many-parameter fits 


e e + , fi fi + , t t + 
e~e + , z < 


95% 
68% 


(0;17.03)M" 1 (0;5.06)M 
(0;27.61)M" 1 (1.68; 7.83)M 


(0;7.40)M 
(0; 12.00)M 


(0; 12.607) 
(0; 33.1288) 



with the SM values of the Z couplings 

SM _ e ( T 3/ ~ 2 Qf Sin2 9 w) SM _ e T 3f 

J sin 6w cos 9w sin 9w cos 9w 

where e is the positron charge, Qf is the fermion charge in the units of e, 
^3/ = 1/2 for the neutrinos and w-type quarks, and T 3 / = —1/2 for the charged 
leptons and <i-type quarks. 

The results of fits of the Z' couplings to the SM leptons obtained from the 
analysis of LEP experiments are adduced in the Tables [5][6l Remind that due to 
the universality of the axial- vector coupling a f the same estimates also hold for 
quarks. First of all, one parameter fits of LEP experiments as well as the many- 
parameter fit for the e + e- backward bins show the hints of the Z' boson at the 
1-2(7 CL. Due to this fact, the fits allow to determine the maximum likelihood 
values of Z' parameters. In spite of uncertainties, these values can be used as a 
guiding line for the estimation of possible Z' effects in the Tevatron and LHC 
experiments. The maximum likelihood values are given in TableJSJ As it is seen, 
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Figure 8: The maximum likelihood values and the confidence intervals for the 
Z-Z' mixing angle (9q) and the axial- vector couplings to the SM fermions (off) 
by the LEP 2 data. The values excluded at 95% CL by the many-parameter fit 
of e + e~ — > are shown in dark gray. The results of fits based on the one- 
parameter observables are shown in light gray for e + e~ — > (i + /i~,t + t~ (left) 
and for e + e~ — > (right). The maximum likelihood value arc plotted as 

the dashed lines. The dotted lines correspond to the maximum likelihood values 
obtained from the LEP1 data. 



different fits obtained for different processes lead to the comparable values of 
the Z' parameters. 

In Table [6] we present the confidence intervals for the fitted parameters. It 
gives a possibility to estimate the uncertainty of the Z' couplings as well as 
the lower bounds on the parameters. The results of fits are also shown in Figs. 
[8][9] To summarize them we note that the data of the LEP2 experiments are 
compatible at 1-2 a level with the existence of the not heavy Z' boson. 

Now we compare the above results with the ones of other fits accounting for 
the Z' presence. As it was mentioned in Introduction, LEP collaborations have 
determined the model dependent low bounds on the Z' mass which vary in a 
wide energy interval dependently on a model. The same has also been done for 
Tevatron experiments. The modern low bound is mz> > 850 GeV. It is also 
well known that though almost all the present day data are described by the 
SM [U H [31 H [SI [55], the overall fit to the standard model is not very good. In 
Ref. [9] it was showed that the large difference in sin 2 0*5"* from the forward- 
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~ee at 95% C.L. : 





,TeV 



Figure 9: The maximum likelihood values and the confidence intervals for the 
vector coupling to electron current (v' e ) by the LEP 2 data. The values excluded 
at 95% CL by the many-parameter fit of e + e~ — > l + l~ arc shown in dark 
gray. The left panel represents the results of fits based on the one-parameter 



observable for e + t 



~e . The right panel shows the la CL area for the 



many-parameter fit of backward bins of e + e 
value are plotted as the dashed lines. 



The maximum likelihood 



backward asymmetry A h ^ b of the bottom quarks and the measurements from 
the SLAC SLD experiment can be explained for physically reasonable Higgs 
boson mass if one allows for one or more extra U (1) fields, that is Z'. A specific 
model to describe Z' physics of interest was proposed which introduces two type 
couplings to the hyper charge Y and to the baryon- minus- lepton number B — L. 
Within this model by using a number of precision measurements from LEP1, 
LEP2, SLD and Tevatron experiments the parameters ay and clb-l of the 
model were fitted. The presence of Z' was not excluded at 68% CL. The value 
of ay was estimated to be of the same order of magnitude as in our analysis and 
is comparable with values of other parameters detected in the LEP experiments. 
The erroneous claim that ay is two order less then the value derived from our 
Table |4] is, probably, a consequence of some missed factors. The upper limit on 
the mass was also obtained mz' < 2.6 TeV . 

These two analyzes are different but complementary. A common feature of 
them is an accounting for the Z' gauge boson as a necessary element of the data 
fits. The results are in accordance at 68-95% CL with the existence of the Z' 
which has a good chance to be discovered at Tevatron and/or LHC. 
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